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1* INTRODUCTION 


With advancements in space flight technology, a new generation of 
aerospace vehicles will soon come into service. Such a vehicle, designed 
to have lifting capability, can operate as a satellite for an extended 
period outside a planetary atmosphere, and upon accomplishing its mission 
can enter the atmosphere and use its aerodynamic maneuverability to reach 
a prescribed region before performing an approach and landing on an air- 
field like an ordinary airplane. 

The portion of the trajectory, starting at the top of the sensible 
atmosphere and ending at a point at low altitude and low speed where ap- 
proach and landing procedures can be initiated — the entry trajectory — is 
the most critical portion of the flight path. It is during this phase of 
flight that the speed is reduced from orbital speed at high altitude to 
subsonic speed at low altitude. During this rapid decrease in the kinetic 
energy, the deceleration, dynamic pressure and heating rate all vary 
markedly. It is important for preliminary design or mission planning to 
have accurate, yet simple formulas expressing the behavior of different 
trajectory variables and physical quantities associated with the trajectory. 

For planar entry trajectories, there exist several analytical theories 
with various degrees of sophistication. Mathematically, one of the best 
theories is Chapman’s theory for analyses of planetary entry CRef. 1). 

In Chapman’s formulation, the pair of equations for planar entry into a 
planetary atmosphere is reduced to a single, ordinary, nonlinear differen- 
tial equation of the second order by disregarding two relatively 
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small terms and by introducing a certain mathematical transformation. The 
equation is integrated numerically for different cases of entry and the 
functions obtained, called the Z-f unctions, are tabulated (Ref. 2). An 
outstanding feature in Chapman’s theory is that the tabulated data of the 
Z-functions are linlversal in the sense that every solution determines the 
motion and heating of a vehicle of arbitrary weight, dimensions, and shape 
entering a planetary atmosphere. 

A major deficiency in Chapman’s theory is that, because of his two 
main assumptions, the equations derived are only approximate, and the ap- 
plications are restricted to entry trajectories with small flight path 
angle or small lift-to-drag ratio. 

Recently, these restrictive assumptions have been successfully re- 
moved. A set of exact equations for three-dimensional entry trajectories 
has been developed using modified Chapman variables to perform the mathe- 
matical transformation (Refs. 3-6). 

In this report the method of directly matched asymptotic expansions 
is applied to the exact equations for three-dimensional entry in terms of 
the modified Chapman variables, resulting in an accurate analytical solu- 
tion. The two-regime approach of directly matched asymptotic expansions 
has proved to be feasible and effective in previous, more restricted, ap- 
plications (Refs. 7 - 10). Now, applied to these exact equations, a power 
ful, useful solution appears. 

The equations of motion for three-dimensional flight about a nonrotat 
ing spherical planet and inside of its atmosphere, assximed to be at rest, 
are derived in Section 2. The exact dimensionless equations using a set 
of modified Chapman variables are derived in Section 3. In Section 4, 
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these equations are transmuted into a form most suitable for an analytical 
integration using the method of matched as)onptotic expansions. The two 
solutions, one valid in the outer region where the gravitational force is 
predominant, and the other valid in the inner region where the aerod 3 mamic 
force is predominant, are obtained. The two solutions are matched directly 
and the composite solution, uniformly valid ever 5 rwhere, is constructed. 

2. BASIC EQUATIONS OF MOTION 

The vehicle is considered as a mass point, with constant mass m , 
moving about a nonrotating spherical planet. The atmosphere surrounding 
the planet is assumed to be at rest and the central gravitational field 
is the usual Inverse square force field. 

r(t) = position vector 

^ ( 2 . 1 ) 

V(t) = velocity vector 

The initial reference frame OXYZ is the planet-fixed system with 
0 at the center of the gravitational field. The OXY plane is referred 
to as the equatorial plane (Fig . 1) . 

The position vector r is defined in this planetocentric system by 
its magnitude r , its longitude 0 , measured from the X-axis, in the 

equatorial plane, positively eastward, and its latitude > measured 
from the equatorial plane, along a meridian, and positively northward. 

The velocity vector V is expressed in terms of its components in 
a rotating coordinate system Oxyz such that the x-axis is along the 
position vector, the y-axis in the equatorial plane positive toward the 
direction of motion and orthogonal to the x-axis, and the z-axis completing 
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Fig. 1. Coordinate Systems. 
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a right handed system (Fig. 1). Let y te the angle between the local 
horizontal plane, that is the plane passing through the vehicle and paral- 
lel to the Oyz plane, and the velocity vector V . The angle y is 
termed the flight path angle and is positive when ^ is above the horizon- 
tal plane. Let ijj be the angle between the local parallel of latitude 

and the projection of V on the horizontal plane. The angle ip is termed 
the heading and is measured positively in the right-handed direction about 
the x-axis. Let i , j^ , and 5 be the unit vectors along the axes of 

the rotating system Oxyz . We have 

? = r i (2.2) 

and 

V = (Vsiny) i + (VcosYCosi|») j + (Vcosysimi)) $ C2.3) 

The system Oxyz is obtained from the system OXYZ by a rotation 0 
about the positive Z-axis, followed by a rotation <|) about the negative 
y-axis. Hence the angular velocity « of the rotating system Oxyz is 

Q = (sinij) i - (^) j + (cos4> k (2.4) 

We deduce the time derivative of the unit vectors i , j , and k 
with respect to the inertial system OXYZ 

f| - M - (cos* J . (fi) J 

^ = - (cosi]) ^) i + (sin(j) ^) k 

| = (f)t- (sin^f)! 


(2.5) 
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The equations of motion of the vehicle are 


and 


dt ■ 


V 


dt 




+ g 


( 2 . 6 ) 


(2.7) 


where g is the acceleration of gravity and the forces L and D are 
the lift and the drag. Expressed in components along the rotating axes, 
g is simply 

g = - g(r) t (2.8) 


The drag D is always opposite to the velocity vector, while the lift L 
is orthogonal to it. Hence, based on Eq. (2.3), we have immediately for 
the vector D 


$ = - (Dsiny) ^ - (DcosycosiJ)) ^ - (DcosysintJ;) ^ (2.9) 


In planar flight, the vector L is in the (r , V) plane and there 
is no lateral aerodynamic force. By control action, if we rotate the 
vector L about the velocity vector V we create a lateral component of 
the lift force that has the effect of changing the orbital plane. To re- 
solve the lift L into components along the rotating axes, we refer to 
Fig. 2. The vertical plane considered is the (r , V) plane. Assume the 
vector L is rotated out of this plane through an angle a . The angle 
a which is the angle between the vector L and the (r , V) plane will 

be referred to as the roll, or the bank, angle. The lift is decomposed 

' > 

into a component Lcosa in the vertical plane and orthogonal to V and 
^ 

a component Lsina orthogonal to the vertical plane. Let x , y , 
and z’ be the axes from the position M of the vehicle, parallel to 
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Fig. 2. Aerodynamic Forces 
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the rotating axes x , y , and z . Let 


be 


the axes from the point M , along the direction of Lcosa , V , 


and Lsina , respectively. The system Mx^y^z^ is deduced from the 
system Mx’y’z’ by a rotation in the horizontal plane, followed by 
a rotation y in the vertical plane. Hence, we have the transformation 
matrix equation 


or 


0 0 
cos4^ -sinijj 


cosy siny 0 
-siny cosy 0 


U'. 

1 

1^0 sin4^ 

cosij; J 0 

0 

ij 


x’ 


cosy 

siny 

0 


1 

y’ 

= 

- sinycosii^ 

cosycos^^ 

~ sin4^ 


^1 

z’ 


- sinysini/^ 

cosiysinilj 

cos^ 

1 

^ / 




( 2 . 10 ) 


Since the components of L in the Mx^y^z^^ system are x^ = Lcosa , 
y^ = 0 , = Lsina , we deduce the components of L along the system 

Mx’y’z’ , or what is the same, along the rotating system Oxyz 


L = (Lcosacosy) i - (LcosasinycosiJj + LsinasiniJ;) j - 

- (Lcosasinysinii^ - Lsinacosij^) k 


( 2 . 11 ) 


Now, if we take the derivative of r , as given by Eq. (2.2), using 
Eq, (2.5) for the derivative of i , we have 

By substituting into Eq. (2.6) and using Eq. (2.3) for the components of 
V , we have three scalar equations 


dr 

dt 


9 

Vsiny 


(2.13) 


Vcosycos^ 

rcos(|) 


(2.14) 


d^ _ Vcosyslnil) 


(2.15) 


These equations are the kinematic equations. 

On the other hand, if we take the derivatives of the velocity vector 
V , as given by Eq. (2.3), using the Eqs. (2.5) for the derivatives of 
the unit vectors i , j , and k , and subsequently the Eqs. (2.14) 
and (2.15) for d6/dt and d(Ji/dt , we have 

dV r . dV . „ dY 2 , ->• 

+ Icosycosij; ^ - Vsinycos4' ^ - Vcosysini)) 


“(siny - cosysin\J;tan(|)) ] j (2.16) 


+ [cosysinij; “ Vsinysini|j 4“ + Vcosycosi|; 

at at at 

2 ^ 

+ — ^^^(sinysini|) + cosycos ij;tan(j))] k 


By substituting into Eq. (2.7), and using the Eqs. (2.8), (2.9), and (2.11) 
for the components of g , D , and L , we have three scalar equations 


. dV . „ dy V^cos^"^ 

siny ^ + Vcosy ^ ^ 


^ cosacosy - “ siny - g (2.17) 


cosy ^ - Vsiny ^ - Vcosytanif^ ^ ^ (siny - cosysinil^tanc})) 


(2.18) 


= - — cosasiny - ^ sinatanii^ - ~ cosy 
m ' m ^ m 

cosf ^ - VslnY ^ 

dt ' dt tanii) dt r ^ ' sinijj 

L , . Lsina D 

= - — cosasiny + — — ; cosy 

m ' mtani^ m 


(2.19) 
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Solving for the derivatives dV/dt , dy/dt , and di|;/dt , we have 


three scalar equations. 

the force equations. 



^ =-^-gsinY 
dt m ® ^ 

(2.20) 

V ^ 
^ dt 

L 

= — cosa - gcosY + — COSY 
m r 

(2.21) 

^ dt 

2 

= - — cosycosij^taniii 

mcosY r 

(2.22) 


The six equations, Eqs. (2.13) - (2.15), and (2.20) - (2.22) are the 
exact equations of motion for flight over a spherical, nonrotating planet 
with its atmosphere at rest. 


3. DIMENSIONLESS EQUATIONS USING MODIFIED CHAPMAN VARIABLES 


In the equations of motion derived in Section 2, the aerodynamic 
lift and drag are now assumed to be 


B - fpsc/ 


(3.1) 


where and are the lift and the drag coefficients, assimied inde- 

pendent of the Mach number and the Reynolds number in the hypervelocity 
regime. The density of the atmosphere, p , is assiamed to be a known 
function of the radial distance r . For each flight program, the con- 
trol functions C , C and a are prescribed functions of the time, 

Li D 

and the integration of the system of six equations of motion requires pre- 
scribing the six initial values for the state variables. In addition, for 
a given vehicle, the parameter S/m must also be given. 
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For a planar trajectory with constant angle of attack, with the ini- 
tial point taken as the origin for measuring the range, that is with 
9^ = 0 , it is seen, by inspection of the equations of motion, that the 

following quantities must be prescribed: 

a/ For the vehicle and flight parameters 

SC-/m and SC_/m C3.2) 

1j D 

or equivalently 

SC^/m and C3.3) 

b/ For the Initial conditions 

r^ , and (3.4) 


By a very ingenious coordinate transformation. Chapman has introduced 
two dimensionless variables defined as (Ref. 1) 


u 


Z 


VcosY 

psc^ 

2m 



(3.5) 


where 3 is the atmospheric density height scale defining the atmosphere 
through the differential relation 

^ = - e(r)dr (3.6) 

P 

Through the transformation (3.5) and through some simplifying assumptions, 
Chapman has reduced the basic equations of motion to a single nonlinear 
differential equation with Z as the dependent variable and u as the 
independent variable. Chapman’s equation is 

— d /dZ 1 - u^ h ^ 3 ^ 

u * 3 ) 33 — cos Y + vpr — cos y = 0 

du du u u Z D 


(3.7) 
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with the flight path angle y obtained from 


^ ^ siny (3.8) 

du u 

The varying quantity r enters the equations as the product 8r . In 

the lower altitudes of planetary atmospheres, where aerodynamic forces 

are effective, the quantity 3r oscillates about a mean value. The 

flight parameter appears in Chapman's equation, Eq. (3.7), as the lift-to- 

drag ratio C^/C^ • Furthermore, for entry trajectories, since the ini- 

1j d 

tial value of Z is nearly zero, only the Initial values u^ and 
need be specified for the integration of Chapman ^s equations. Hence, for 
a given atmosphere, with 3r considered as constant (e.g., for the Earth 
3r - 900) , Chapman can integrate numerically his equations and tabulate 
the results for each set of values of • These 

Tables of the Z functions (Ref. 2) can be used for any vehicle of arbi-- 
trary weight, dimensions and shape entering the specified planetary atmo- 
sphere. Furthermore, all the physical quantities during entry, such as 
deceleration, dynamic pressure, heating rate, and heat transfer rate, can 
be easily obtained as simple functions of the variables Z , u and y 
The major deficiency of Chapman theory is that, because of his two 
main assumptions, namely that 

a/ the percentage change in distance from the planet center is 
small compared to the percentage change in the horizontal component of the 
velocity, that is 


dr 

r 


d(Vcosy) 

Vcosy 


(3.9) 


b/ the quantity tany is small, that is 


and 
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L ^ 

— tany 


« 1 


(3.10) 


the validity of Chapman numerical analysis is restricted to trajectories 
with small flight path angle, or trajectories with small lif t-to-drag 
ratio. 

These restrictions have been removed by deriving the exact equations 
for three-dimensional reentry using a set of modified Chapman variables 
(Refs. 3-6). Chapman’s assumptions (a) and (b) , Eqs. (3.9) and (3.10), 
imply that the independent variable u is monotonically decreasing. This 
is only true for the last portion of the trajectory since at high altitude 
u is oscillatory, and in the limit, for flight in the vacuum, u is 
purely periodic (Ref. 5). To avoid this difficulty the following dimension- 
less variable is introduced as the independent variable 


s 


^ V 

/ 3 COSY dt 
0 ^ 


( 3 . 11 ) 


The variable is strictly increasing as long as cosy >0 , a condition 

which is always satisfied for entry at constant lift-to-drag ratio, and 

for all physically realistic entry trajectories. For the equations, 

2 2 
V cos y 
u = ^ 

(3.12) 



Expressed in terms of the original Chapman variables, we simply have 

Z = £ , u = (3.13) 

u 

The new variables lead to a set of differential equations in a simpler 
form allowing a complete qualitative discussion of the three-dimensional 
reentry trajectory. 
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We shall present below a slightly different derivation of the exact 
dimensionless equations from the one given in Ref. 5. 

First, by eliminating the time, by dividing Eqs. C2.14), (2.15), and 
(2.20) - (2.22) by Eq. (2.13), we have 

d6 _ cosil' 
dr rcos^tany 


dj) _ sinij) 
dr ” rtany 


dr msiny 


^ = 
dr 


pSC^cosa 
2ms iny 


dij; 

dr 


pSC sina 
L 

2msinycosy 


2g 


+ -i- 

V^tany 

cos4^tan(j) 

rtany 


Next, from Eq. (3.11), using Eq. (2.13), we have 


ds _ ^ _ 1 

dr dt dr rtany 


(3.14) 


(3.15) 


Hence, using this equation, we can rewrite the Eqs. (3.14) with s as 
independent variable 


d6 _ cos4/ 
ds cos(j) 


^ = sin^ 



^ = 
ds 


di|j 

ds 


rpSCjjV^ 

mcosy 


- 2grtany 


rpSC- cosa 

L + n _ EL) 

2mcosy '' ^ 2 ^ 

rpSC sina 

Li 


(3.16) 


2mcos y 


- cos4^tan(j) 
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Using the Chapman formulation, the variable V is replaced by the var- 
iable u , while the radial distance r is replaced by the variable Z 
By taking the derivative of u , as defined by the first equation (3.12), 
with respect to s , using an inverse-square law for the acceleration of 
the gravity, we have 


2 2 2 2 2 
du _ cos Y dV _ 2V sinycosY dy ^ V cos y dr 

ds gr ds gr ds 2 ds 

8 ^ 


(3.17) 


Using the definition (3.12) for u and Z , with the appropriate deriva- 
tives from Eqs. (3.15) and (3.16), we have the differential equation for 


du ZUr- . / Lv . . siny" , 

XT ■[! + (tt-) cosotany + 


ds 


cosy 


2 /^ Z 


(3.18) 


On the other hand, by taking the derivative of Z , as defined by the 
second Eq. (3.12), with respect to s , we have 


ds 


= - Br 


'■ p3 dr 


2Br 



) Ztany 


(3.19) 


where B' = dB/dr 

As discussed in Ref. 3, the term in parentheses in equation (3.19) is 

1 3 * 

for the locally exponential atmosphere of equation (3.6), 1 - — T 

2B 

2 

For an isothermal atmosphere Br is constant and this term becomes 


1 - 


1 dp 3 

pB dr 2Br 

1 


. For a strictly exponential atmosphere this term is 


. Finally, for an atmosphere with Br constant, this term is 


1 dp 


. In all these cases the direct dependency on r is of 


pB dr Br 

the form Br . In fact, for Earth, where Br is about 900 , the term 


in question will be unity to a high degree of accuracy. As analyzed in 
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Refs. 3 and 5, the equations with gr as a constant may still justifiably 
be considered exact. 

Finally, the differential equations for y ^ > written in terms 

of the dimensionless variables are, respectively 


and 


ix == Z j. A) 

ds cosy^^'C^ 


COSO + -2221(1 

/6r Z 




Q 2 

di|) _ /g7 g^na cos ycosij;tan(|) j 

cos^y /gr Z 


(3.20) 


(3.21) 


In summary, the exact equations of motion for three dimensional flight 
in an isothermal atmosphere, using modified Chapman variables are 


ds 


- 3r(l - -^) Ztany 


^ . . ML?S(i + itanY + 

ds cosy 2 /^ 2 

iX = hx + -^°-g-X(i - co^)3 

ds cosy 2 "" 

d9 _ cosijj 
ds cos(j) 


(3.22) 


ds 


= sinijj 


where 


dijj _ /gr Z .Q _ cos ycosi|)tan(j) ^ 

ds 2 „ 

cos y vgr Z 


X = — COSO 


6 = — sina 


(3.23) 


It may be observed that, although the equations are derived for three- 
dimensional flight, the first three equations are decoupled from the last 
three equations. Hence, by using only the assumption of constant value 
for gr , for constant lift-to-drag ratio and constant bank angle, X 
is the only flight parameter that needs to be specified. For each X , 
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with specified initial conditions ~ 0 » the first 

three equations (3.22) can be Integrated and Tables of the Z functions 
can be generated. Since Chapman’s two restrictive assumptions have been 
removed, the new Tables are valid for three-dimensional flight without 
any restriction concerning the flight path angle and the lift-to-drag 
ratio. 

Reference 5 presents a qualitative and quantitative analysis of these 
equations. It is the purpose of the present work to obtain an analytic 
solution to the complete set of the exact equations of motion. These 
equations are integrated by the method of matched asymptotic expansions. 
The results of investigation corroborate the assessment in Ref. 5 that 
the assumption of a constant 6r is very accurate. Hence, the equations 
(3.22) with gr = constant can be considered as the exact equations for 
reentry, and Tables of the Z functions based on the new equations should 
provide accurate data for analyses of planetary entry. 


4. SOLUTIONS BY DIRECTLY MATCHED ASYMPTOTIC EXPANSIONS 

In this application of the method of matched as 3 miptotic expansions, 
the solutions are obtained separately for an outer region, near the 
vacuum, where the gravity force is predominant, and for an inner region, 
near the planetary surface, where the aerodynamic force is predominant. 
Hence, the altitude is the appropriate independent variable selected for 
the integration. 

Let y be the altitude and let subscript t denote the reference 
altitude, taken at sea level. Then 

r = r + y = r (1 + h) (4.1) 

t t 
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where the dimensionless altitude h is defined as 


h = ^ 


The differential relation between s and h is 


(4.2) 


ds = 


dh 


Cl + h) tany 


(4.3) 


For the integration, we adopt a strictly exponential atmosphere, but the 

general method can be applied to any more realistic atmosphere such as, 

for example, the one proposed in Ref. 7. For an exponential atmosphere 

h 


-3y 

P = P.|. e e 


where 


e = 


art 


(4.4) 


(4.5) 


Since the constant 3r^ is large, e.g., for the Earth atmosphere 
3r^ = 900 , the parameter e is a small quantity. By the definition 

(3.12) of Z 


2 . /UThT j E 

2 ra 3 e ® 


(4.6) 


We define the ballistic coefficient 


SC 


B = 


D^t 


2m3 


(4.7) 


For each vehicle, B is specified and the variable Z is obtained from 

Z = e (4.8) 

By this relation, the first Eq. (3*22) can be deleted, and we write the 
other equations with the dimensionless altitude h as independent vari- 


able 
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_ h 

du u 2Bu(l 4- Xtany) e 

dh Cl + h) esiny ° 


ia = _ 

dh 


— 9 — (1 - _ 

Cl + h) u '' 


AB 


h 

e 


d6 _ cosij) 

dh (1 + h) cos(|itanY 


(4.9) 


dj> sinij; 

dh Cl + h) tany 


h 

0 

dip _ cQsil^tanc}) ^ B6e 

dh (1 + h) tany esinycosy 


We have defined 


q = cosy (4.10) 

Also, we recall the definition of the flight parameters 

X = — cosa , 6 = — sina (4.11) 

The Eqs. (4.9) are in a suitable form for numerical integration 
for flight inside an atmosphere. For an analytical solution of the re- 
entry trajectory using the method of matched asymptotic expansions, we 
shall use a more convenient form using some elements of the orbit as intro- 
duced in celestial mechanics, since these elements are constants of the 
motion for flight in a vacuum. 

As seen in Fig. 3, if I is the inclination of the plane of the os- 
culating orbit, that is the (r , V) plane, Q the longitude of the as- 
cending node, and a the angle between the line of the ascending node 
and the position vector, we have the following pertinent relations from 
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spherical trigonometry 


cos((>cosi}; = cosi 

" ISl 

cosa = cos(J)cos(0 - Q) 

These relations are independent. We can easily deduce 


sincj) 

sinij; 

sinij^ 


sinisina 

tan<{) 

tana 

sinIcosCG “ 


(4.13) 


Using these relations, we replace the variables 0 , ^ , and i|; by 

the new variables a , , and I . The Eqs. C4.9) now become 

_ h 

du u 2Bu(l + Xtany) e 

dh " (1 + h) esiny ® 


2 - — 
ia g n _ i-') _ M e ^ 

dh (1 + h) u e ® 

tl 

da 1 BSsina e 

dh ~ (1 + h) tany etanisinycosy 

_ h 

df^ B6sina £ 

dh ~ esinisinycosy 

_ h 

dl B6cosa e 

dh esinycosy 


(4.14) 


The Eqs. (4.14) are most suitable for an integration using the method 
of matched asymptotic expansions. We notice that, once the elements a , 

, and I are known, we obtain the original variables 9 , <j) , and \p 


from 
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tan(0 - f2) =« cositana 

sin<}) = slnlsina (4.15) 

tanij^ = cosatani 


4.1 Outer Expansions (Keplerian Region) 


The Eqs. (4.14) are expressed in terms of the outer variables. The 
outer expansions are introduced to study the limiting condition of the 
solution in the region near the vacuum where the gravitational force is 
predominant. They are obtained by repeated application of the outer limit, 
which is defined as the limit when e 0 with the variable h and 
other dimensionless quantities held fixed. 

We assume the following expansions 


u = UQ(h) + £u^(h) + . . . 

q = <1 q(Ii) + £q-j^(h) + . . . 

a = eUj^(h) + . . . 

Q = fJQ(h) + eS^^(h) + . . . 

I = lQ(h) + el^(h) + . . . 


(4.16) 


By substituting into Eqs. (4.14) and equating coefficients of like power 
in e , the differential equations with zero order of e are 

du^ u.^ 


dh 

dq^. 

dh" 

da 


(1 + h) 


-^0 g _ :i0 

(1 + h) Uq^ 


0 


dh (1 + h) tany. 


(4.17) 


diJj. 

dtT 

dh 


= 0 


= 0 



The solution of this system is 
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Uq( 1 + h) = 


q,^ 1 


u. 


1 + 


Cj^C2 cosCOq - C^) 


(4.18) 


% ^4 


^0 = S 


where the are constants of integration. The first and higher order 

solutions are all equal to zero because at high altitude, in the limit 
the atmospheric density is zero and the motion is Kepler ian. 


4.2 Inner Expansions CAerodynamic-Predominated Region) 

The inner expansions are introduced to study the limiting condition 
of the solution near the planetary surface where the aerodynamic force is 
predominant- They are obtained by repeated application of the inner limit, 
which is defined as the limit when e 0 with the new stretched alti- 
tude 

h = - (4.19) 

e 

and the other dimensionless quantities held fixed. 

We assume the following expansions 

u = u^Ch) + eu^(h) + . . . 

q = ^0^^^ ^ + . . . 

a = otgCh) + ea^(h) + . . . (4.20) 

Q = f^gCh) + e^^^(h) + . . . 

I = IpCh) + elj^Ch) + . . . 
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By substituting into Eqs. (4.14) and equating coefficients of like powers 
in e , the differential equations with zero order of £ are 




dh 

dS^ 

dh 

dl 

dh 


0 

0 


2BUq( 1 + AtanYQ> 

e 

siny 

, _ -h 
- XBe 


B5sina, 


tanlQsinYQCOSYQ 
BdsinuQ 

siniQsinYQCOSYQ 
B6cosa, 


-h 


0 


-h 


sxnYQCosYQ 


(4.21) 


The solution of this system is 


~ ~2 ^^0 
Uq = Cj_qQ exp[- — ] 


q^ = XBe ^ + C, 


sinaQsinlQ = sinC^ 


(4.22) 


cosa^ = 


cosC^cosCC^ - ^q) 


n . "^ 0 . 


cosIq = cosC 2 Cos{y log! tan (^ + + C^} 


where the C are constants of integration. 
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4.3 Asymptotic Matching and Composite Expansions 

The constants of integration Cj in the inner expansions will be 
determined by matching with the outer expansions. In this problem, 
matching is accomplished by expanding the inner solutions for large h , 
expressing the results in terms of the outer variables and matching with 
the outer solutions for small h 

The outer solutions, Eqs. (4.18), become for small h 


‘lo 

“o 


2 - 


-1 ^ 

cos ( — • = : ) + C, 


(4.23) 



On the other hand, the inner solutions, Eqs. (4.22), become for large h 

~ ~ ~2 2 - 1 ~ 

Uq = C^C 2 exp[- Y cos 

«0 ' =2 

sinagSinlp = sinC^ (4.24) 

cosUq = cosC^cosCC^ - ^^q) 

cosIq = cosC^cos{y log [tan (^ + ^ ^^ 2 )] + C^) 
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Matching Eqs. (4.24) with Eqs. (4.23) provides the constants In terms 

of the constants . We have 


r2 -1 


= C2 - sxply cos / 2 - ] 


/ 

\ - 2 


-1 '^1 ■ ^ 

sin Cg = sinC^sln[cos ( ^ — — ) + C.,] 


A 


A^2 


(4.25) 


-1 -1 ^ 

C, = cos {cosicos ( ) + C^]/cosC^} + C, 

4 ) ^ J 4 

A - 


/ ^ 


= cos 


^[cosC^/cosC^] - y log[tan(^ + -j cos ^ / 2 -* ^^^2 


Hence, the constants are expressed explicitly in terms of the con- 

stants C . Substitution into Eqs. (4.22) gives the inner solutions. 
J 

It is convenient to use the following notation to write these solutions 
in a symmetric form. Let 




cos 


. / 2 - CjCj 


-1 ^ 

sin ()>^ = sinC^sin[cos ( — ) + C.^] 


crc 


(4.26) 




-1 -1 A ■ ^ 

0. = cos {cosicos ( ) + C^] /cOS(t)j.} + C, 

* / — 3*4 

/i - 


I* “ s 
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The constants with subscript * are explicit functions of the constants 
Cj . Then the inner solutions are 

— = - Yq)] 

* cos Y* 

cosYq = cosy* + XBe ^ 

sina^sinlQ = sin^)^ (4.27) 


cosOq = cos(f)^cos(6^ - 

cosIq cosl^ 

7— ) “ cos ( T“) 

I COS(J)^' COS(|)^ 

j 

I 

i 


- 1 . 
COS C 


J log [ tan (|- + ^)/tan(|- + -^)] 


The Eqs. (4.27) show that during the phase of aerodynamic force 
predominant turning, the latitude (j) and the longitude 0 remain constant. 
The last equation gives the change in the heading ip during that phase. 

The composite expansions, valid ever}rwhere, can be constructed by 
the method of additive composition. The additive composition is 
obtained by taking the sum of the inner and the outer expansions, 

Eqs. (4.27) and 4.18), and subtracting the part they have in common 
(the inner limit of the outer expansions or the outer limit of the 
inner expansions), Eqs. (4.23) or (4.24). Thus, for the variables 
u and Y> using subscript c for the composite solution. 


u 




, cos Yr. n ^ 

^ exp[f(Y* - y^)] 


(1 + h) ■ 2 ^0 

cos Y* 


(4.28) 


cosy = cosy^ / r » r- 

J 2cos Y^(l + h) + (u^ - 2cos Y*) (1 + 


+ XBe 


h 

e 
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For the angular variables a , Q and I , we have 

-1 S ’ 1 

a = - cos ( ) 

c 0 U J j ^ — 


^ - C, 

c 0 0 4 


(4.29) 


= ^0 + "O - S 


Hence, from equation (4.23), we have immediately , and 




—1 6 n ^0 n 

cosi = cos(() . cos{cos ( — ) + -:r log I tan (y + — )/tan(y + y)]} 

C ^ COSq)^ A 4 Z 4 z 


(4.30) 


For the angle , we can use the second Eq. (4.12) to have 


-1 


tanij)^ 


^ = 0* “ sin [- — y] 

c * tani 

c 


(4.31) 


where I is given by Eq. (4.30) 


Finally, the angle is given by 

sin(J)^ cosy^ 

a = sin [ . y " J + cos [ ' - ■ - -y 

c sxnt / ty 2 

^ /uf + (1 - 2u.) cos Y 


1 + h 
* 


- 1 )] 


(4.32) 


_1 cosY*(u* - 1) 

- cos [- 


+ (1 - 2u.) cos^Yj 


The composite solutions are expressed explicitly in terms of the 

constants of integration » Y* » <!>* » 0* and . For the 

computation in terms of the Independent variable h , the angle Yq 

is first calculated from the second Eq. (4.27). Subsequently, we have 

u , Y and I , and then and a 

c c c c c 
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4.4 The Composite Solutions In Terms of the Initial Conditions 

For the initial conditions to be satisfied identically, the five con- 
stants of integration Cj , or equivalently the five constants with sub- 
script * , as defined by Eqs. (4.26), are to be evaluated by using the 

composite solutions. Let the conditions at h. be 

u=u. ,Y=Y. ,a=a. ^ Q = Q. ,1=1. (4.33) 

x x’ X X X 

Using these conditions in the composite solutions, the constants u^ , 

Y* > > and are obtained upon solving a set of transcen- 

dental equations which can only be done ntjmerically. Another obstacle 
arises when, as is a common practice, in order to reduce the number of 
prescribed initial values, one takes the initial (r^ , V^) plane as the 
reference OXY plane. In doing so we have = 0 , and = 0 , 

but when 1=0 , the longitude of the ascending node is not defined 

as evidenced by Eq. (4.31). This singularity can always be avoided by 
rotating the OXY plane through a fixed and arbitrary angle, say 45° 
about the r, axis. Then the initial condition at h. is 

X 1 

u = , Y = Yi , a. = 0 , 0. = 0 , = 45° (4.34) 

The equivalent condition for the variables 6 , c(i and >p is 

= 0 , <^^ = 0 , = ‘is- (4.35) 

Finally, it should be noted that the composite solutions for u and 
Y , Eqs. (4.28) remain unchanged for the planar case. For the planar 
case, with the motion along the equatorial plane, the variable a is the 
same as the longitude 6 . The composite solution for 0 can be seen 


to be 
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where is a constant of integration. The three constants of integra- 
tion , Y* and in Eqs. (4.28) and (4.36) for the planar case 

are evaluated using the initial conditions > and 0^ 


5. CONCLUSIONS 

For a spherical, nonrotating planet with a spherically symmetric, 
but otherwise arbitrary, atmosphere the exact equations for three-dimen- 
sional, aerodynamically affected flight have been derived. The equations 
are transformed using modified Chapman variables into a set suitable for 
analytic integration using as 3 nnptotlc expansions. 

The solutions by asymptotic expansion are obtained as an inner ex- 
pansion in the aerodynamically dominated region, and an outer expansion 
in the gravitationally dominated region. 

The inner and outer expansions are matched directly, and a composite 
solution, valid everywhere, is constructed by additive composition. 

This method of directly matched asymptotic expansions has provided 
highly accurate and useful solutions to less general atmospheric trajec- 
tory equations. Ref. 7. In that work extensive numerical calculations 
demonstrated the accuracy compared with exact numerical solutions. 

The behavior of the exact entry equations was examined in Ref. 5. 

The method of directly matching the inner and outer expansions for at- 
mospheric entry trajectories was proven valid for some restricted prob- 


lems in Refs. 8-10. 
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This report has accomplished the wedding of the exact atmospheric 
trajectory equations, using the modified Chapman variables, with the 
method of directly matched asjnnptotic expansions and provides an analy- 
tical solution which should prove to be a powerful tool for aerodynamic 
orbit calculations. 
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